Astronomy & Astrophysics manuscript no. paper 


©ESO 2013 


March 3, 2013 





Rapid growth of gas-giant cores by pebble accretion 



Michiel Lambrechts 1 and Anders Johansen 1 



'Lund Observatory, Department of Astronomy and Theoretical Physics, Lund University, Box 22100, Lund, Sweden 
e-mail: michiel@astro.lu.se 



Received — ; accepted - 



ABSTRACT 



i— H 

o 

CO 



Oh 

Cu 
O 

<N 
> 
O 

m 
o 
co 

in 
o 



The observed lifetimes of gaseous protoplanetary discs place strong constraints on gas and ice giant formation in the core accretion 
scenario. The approximately 10-Earth-mass solid core responsible for the attraction of the gaseous envelope has to form before 
gas dissipation in the protoplanetary disc is completed within 1-10 million years. Building up the core by collisions between km- 
sized planetesimals fails to meet this time-scale constraint, especially at wide stellar separations. Nonetheless, gas-giant planets are 
detected by direct imaging at wide orbital distances. In this paper, we numerically study the growth of cores by the accretion of cm- 
sized pebbles loosely coupled to the gas. We measure the accretion rate onto seed masses ranging from a large planetesimal to a fully 
grown 10-Earth-mass core and test different particle sizes. The numerical results are in good agreement with our analytic expressions, 
indicating the existence of two accretion regimes, one set by the azimuthal and radial particle drift for the lower seed masses and the 
other, for higher masses, by the velocity at the edge of the Hill sphere. In the former, the optimally accreted particle size increases 
with core mass, while in the latter the optimal size is centimeters, independent of core mass. We discuss the implications for rapid core 
growth of gas-giant and ice-giant cores. We conclude that pebble accretion can resolve the long-standing core accretion time- scale 
conflict. This requires a near-unity dust-to-gas ratio in the midplane, particle growth to mm and cm and the formation of massive 
planetesimals or low radial pressure support. The core growth time-scale is shortened by a factor 30-1,000 at 5 AU and by a factor 
100-10,000 at 50 AU, compared to the gravitationally focused accretion of, respectively, low- scale-height planetesimal fragments or 
standard km- sized planetesimals. 
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1. Introduction 

The gas giants (Jupiter and Saturn) and ice giants (Uranus and 
Neptune) in our solar system consist of a dense rocky and/or icy 
core surroun ded by a vary ing degree of hydrogen and helium 
atmosphere dGuillotl l2QQ5b . The strong positive co rrelation be- 
tween stellar metallicity and e xoplanet occurrence (ISantos et al. 
l2004t iFischer & Valentill2005l) is also accompanied by a corre- 
lation between stellar metallicity and the amoun t of heavy ele- 
ments present in the e xoplanetary interior (Guillot et al. 2006; 
iMiller & Fortneyll201 ll) . for objects in the gas giant mass range 
between about 0.3 and 10 Jupiter masses (Mj). Additionally, a 
careful statistical inspection of the planet candidates from the 
Kepler transit survey reveals the evaporation and sublimation 
of the smaller ice and gas giant pl anets to their naked cores 
as they get close to their host star (lYoudinll201lb . After their 
formation, the migration of these massive planets in the later 
stage s of the protoplan etary disc shapes the final planetary sys- 
tem ( Walsh et al. 20lTh. However, reconstructing how ice and 
gas giants form in the first place has proven to be challenging. 

In the disc instability scenario, gravitational instabilities in 
the protoplanetary disc excite d ense spira l arms which fragment 
directly into gas giant planets (lBosslll997l) . The core accretion 
(or nucleated instability) scenario requires the formation of a 10- 
Earth-mass (M e ) solid core, capable of holding on to a gaseous 
atmosphere. When the envelope reaches a mass comparable to 
the core mass, a run-away accretion of the surrounding gas is 
triggered dMizunoll98fl iPollack etalJI 1996b . 

Using the solar system as a template for the end result of 
planet formation is challenged by direct imaging of planetary 
companions to A-stars at wide stellar separations. The system 



HR 8799, for example, contains at least 4 planets separated from 
their host star by 14.5, 2 4, 38 and 68 AU, confined by an in- 
ner and outer debris disc (Maroi s et al . 2010). Best estimates of 
the planetary masses place them all in the gas-giant range. The 
presence of the debris discs reveals that growth to planetesimals 
occurs at wide orbital distances as well. Anoth er example of a 
direct ly imaged gas-giant planet, (3 Pictoris b (lLagrange et al.l 
120101) . orbits the host star at appro ximately 10 AU. F omalhaut 
b detected in reflected visble light ( Kalas et ID [2008), with an 
upper mass below 1 Mj ( Janson et al.l l2012), is located far from 
the central start at approximately 120 AU. LkCa 15b is a newly 
discovered gas-giant planet of about 6 Mj, likely caught in the 
epoch of formation, orbiting at approximately 20 AU around 
a young solar-like star , with an estimated age of only 2 Myr 
dKraus & Irelandll2012h . 

Formation of gas giants by direct gravitational collapse has 
bee n shown to be proble matic. At large distances from the host 
star, iKratter etall (1201 Ol) point out that it becomes increasingly 
difficult to clump gas with masses below the deuterium burn- 
ing limit. Additionally, at smaller stellar separations, gas cools 
too slowly for the spiral arms to fragm ent into bound clumps 
dMatzner & Levidl2005URafikovll2005l) . 

On the other hand, gas-giant formation by core accre- 
tion suffers fro m exceedingly long time- scales at wide stel - 
lar separations (iDodson-Robinson et al.l 120091; iRafikovl 1201 ll) . 
Observations of dust infra-red emi ssion (lHaisch et al.l |200l| ; 
llC urrie et al.l 2009) and disc accretion ( Javawardhana et al.ll 2006) 
limit the lifetime of the gaseous component of the protoplan- 
etary disc to 10 6 - 7 yr. Classical core formation by runaway 
planetesimal accretion is believed to take more than 10 7 yr be- 
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yond 5 AU, where th e planetesimal number densities are low 
(Goldrei ch et al.ll20 04). Planetesimals (> km) get gravitationally 
focused on to the core, but this effect can be significantly reduced 
when scattering events drive up the random velocity component 
of the planetesimals. 

The formation of planetesimals, larger- than-km- sized solid 
bodies bound by self-gravity, is problematic in its own way. 
While classically considered to be the building blocks of both 
rocky planets and gas-giant cores, the formation of solids 
this size remains difficult to explain both theoretically and 
experimentally. Part icle growth beyond cm- sizes by coagu- 
lation is ineffici ent felum & WurrnT 120081: iBrauer et al.l 120081: 
IWindmark et al.l 1201 2l) and radial drift time- scales for m- 
- sized boulders are as short as a hundred orbital time- scales 
(lAdachi et alJl!976b IWeidenschilling 1977). On the other hand, 
one can circumvent this so-called meter barrier with tur- 
bulence induced by th e magnetorotational instability (MRI, 
Balbus & Hawlev 1991), which excites local pressur e bumps, 
ideal regions for dust partic le trapping and growth (Whipple 
119721: iJohansen et"aD l2009al) . In dead zones where the MRI 
does not operate, streaming instabilities ca n destabilize the rel- 
ative motion between gas and particles (|Youdin & Goodmanl 
120051: IJohansen & Youdinl 120071: iBai & Stond 12010b and lead 
to the formation of dense filaments. When the particle den- 
sity becomes sufficiently high, large Ceres - sized planetesimal s 
are formed through gravitational collapse (iJohansen et al.ll2007l) . 
The streamin g instability benefits str o ngly from in c reased disc 
metallicities (IJohansen etaD l2009bt IBai & Stond l2010h . ex- 
plaining partly the higher occurrence rate of exoplanets around 
higher metallicity stars. 

Instead of building up cores of ice and gas giants with plan- 
etesimals, we investigate in this paper the accretion of smaller 
particles, coupled to the gas on approximately orbital time- 
scales. Dust continuum observations of young circumstellar 
discs around low-mass pre-main- sequence stars show growth 
of the domin ant particle size to mm and cm sizes within less 
than 1 Myr (iTesti et al.ll2QQ3t IWilner et al.l2005t iRodmann etafl 
2006). The dynamics of these small particles is influenced by the 
prese nce of the surrounding p rotoplanetary gas (IWeidenschillingl 
119771) . through Epstein drag (lEpsteinll 19241) . 

While drag helps re ducing the ran dom velocities of large 
planetesimals (> 1 km), R afiko v (2004) carefully investigated 
analytically the effect of drag on smaller bodies (< 1 km) 
as assumed products of a collisional cascade. However, he 
did not consider particles coupled to the gas on shorter-than- 
orbital time- scales, excluding the pebble- sized objects seen in 
T Tauri discs. He finds that nearly all fragments settle to the 
midplane of the nebula and that gas drag is efficient enough 
to prevent dynamical excitation, making core formation possi- 
ble within near ly 10 6 yr, as was later c onfirmed in coagula- 
tion models by lKenyon & Bromley! d2QQ9h . Accretion of smaller, 
pebble-sized particle s onto protoplanets was first investigated by 
Joh ansen & Lacerdal (l2010h . who numerically found that peb- 
bles are accreted from the entire Hill sphere, the region roughly 
corresponding to the maximal gravitational reach of the core. 
They identify a prograde particle disc, which could explain 
the spin periods of asteroids and preferential prograde spin of 
large asteroids. The influence of gas drag on the interaction 
of single small bodies and low-mass planets was explored by 
lOrmel & Klahrl (120101) . Analytically, they calculated that proto- 
planets starting from ~10 3 km can efficiently accrete ~cm- sized 
particles wi th impact parameters compa rable to the radius of the 
Hill sphere. lOrmel & Kobayashil 12012) further investigated the 
protoplanet growth stage with a thorough toy model including 



fragments, planetesimals and embryos and stressed the impor- 
tance of the gas disc properties , such as a reduced local headw ind 
and turbulence for fast growth. iPerets & Murray-Clayl(l201ll) an- 
alytically investigated the coalescence of binary planetesimals 
due to drag forces and c ommented on the p ossibility of growth 
through this mechanism. iLyra et a l. (2008a) had already earlier 
ran full disc models of pressure bumps formed near the edges 
of the deadzone. After merely 200 orbits, they observed bound 
embryos with masses similar to the planet Mars, consisting of 
pebble- sized particles. 

In this paper, we investigate core growth from a seed mass 
by gas-drag-aided capture of cm-sized pebbles. In Section 2, we 
describe the physics included in the shearing coordinate frame 
used to numerically model the growth of the core. In Section 3 
we present the results from our simulations and analyse the ac- 
cretion rates for various core masses. We compare our results to 
analytic expressions capturing the essential physics underlying 
the phenomena at hand, namely the sub-Keplerian gas velocity, 
the particle size, the Keplerian shear and the gravitational pull 
from the seed core. The effect of local changes in the pressure 
gradient are analysed and we present the effect of including the 
backreaction of the particles on the gas flow. By extrapolating 
the measured accretion rates, we discuss the formation of gas 
and ice giant cores and derive a characteristic time- scale for core 
formation by pebble accretion in Section 4. We discuss the ap- 
proximations made in this paper and the limitations of our model 
in Section 5. Finally, in Section 6, we conclude that pebble ac- 
cretion can explain rapid gas and ice giant formation in the core 
accretion scenario, even at wide stellar separations. 

2. Physical model 

The growth of a gas giant's core occurs in a protoplanetary disc, 
a gaseous disc in the process of accreting onto the young star. 
Based on the mass distribution in the solar system and assuming 
a mean gas-to-dust ratio or metallicity of 

Z=^=0.01, (1) 

with £ p and £ d enoting the solid (dust+ice) and gas column den- 
sities, [Hayashj| ([19811) constructed the minimum mass solar neb- 
ula (MMSN). He found the radial dependence of the gas column 
density to be 



S = 1700 



-3/2 



gem 



(2) 



with the orbital radius r within 0.35-36 AU. The thin disc is char- 
acterized by the ratio of the gas scale height H to the orbital 
distance 

He. I r X 1 / 4 

« = £L* o.033 (-U , (3) 
r vk VAU/ 

with c s the sound speed of the gas and vk the Keplerian velocity, 
(GM\ l/2 

rQ K = ( — ) . (4) 



v K 



\ r I 



Here Qk is the Keplerian frequency. 

Solids with radii smaller than the local mean free path of 
the gas, R < (9/ 4) A, a re in the Epstein regime of gas-particle 
coupling (lEpsteinll 19241) . They react on a friction time-scale tf to 
changes between the relative particle velocity v and the local gas 
velocity w, 



Vdrag = --{V-U) = " (V-U), 

tf p.R 



(5) 
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Fig. 1. The physical response of a particle reacting to the gas 
flow is set by the friction time tf. A particle of given size R has a 
dimensionless friction time Tf = Q^tf that depends on the orbital 
distance r. The red dash-dotted line marks the distance at which 
the particle size equals 9/4 of the mean free path A of molecular 
hydrogen in the MMSN. Particles with R > (9/ 4) A are located 
in the dashed region, and experience Stokes drag as opposed to 
Epstein drag. For the curves in the Stokes regime, we have ig- 
nored the transition into the non-linear Stokes regime, applicable 
for large particles close to the host star. 
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Fig. 2. Deviation A of the orbital velocity of a gas element with 
respect to an object orbiting with the full Keplerian frequency, 
normalized by the local sound speed, is plotted as function of 
orbital radius r in AU. The bold black line represents the tradi- 
tional MMSN scaling, while the b old red line corresponds to the 
adapted MMSN as presented in IChiang & Youdinl OOlOh . The 
shaded area connecting to thin curves indicates the effect of a 
strong pressure bump of strength SA = -0.04. The adopted stan- 
dard value of A = 0.05 in this paper, is accurate in a region 
around 5 AU, even without a strong pressure bump. 



where R and p. are the radius and material density of the particle, 
while p is the local gas density. For particles in the vicinity of 
the midplane, with z < one can assume pH « Z/ V27T, so 
that the particle size R in the MMSN can be recovered from its 
dimensionless friction time 



Tf = Q K t f 

(also known as the Stokes number) as 



R = 60 cm Tf 



2gcm~ 



-3/2 



(6) 



(7) 



Figure \T\ shows the relation between the orbital radius and the 
particle radius for different dimensionless friction times. Around 
10 AU, a dimensionless friction time of Tf = 0.1 corresponds to 
cm-sized particles, which we will refer to as pebbles. Close to 
the star, the gas density increases sufficiently for the particles to 
enter the Stokes drag regime, where t^ = (4/9)(R/A)Tf scales as 
ocr 5/4 . For a more comp let e description of different drag regimes, 
see e.g. Rafikov (2005) or Youdin (2008). 

The gas component of the protoplanetary disc moves with a 
sub-Keplerian mean velocity, since the force due to the the solar 
gravity is reduced by the radially outwards pointing gas pressure 
force. The azimuthal velocity difference A = Au^jc^ between 
the mean gas flow and a pure Keplerian orbit is given by 



vk 



1 c s dln(P) 
~2v K <9m(r)' 



(8) 



where P = pc 2 s is th e gas pressure and rj i s a measure of the 
gas pressure support (Nakagaw a et al . 1986). In the MMSN, A 



has a weak radial dependency, A « 0.05 (r/AU) 1/4 , as can be 
seen in Figure [2] However, comparison of the MMSN model 
with observed protostellar accretion discs, (e.g.lBell et al. 1 19971) 



(e^Bell 
(Lodders 



and studies of solar nebula metallicities (ILoddersI 120031) have 
prompted updated MMSN mode ls, with a less steep p ressure 
gradient, A = 0.036 (r/AU) 2/7 (IChiang & Youdinl 1201 Oh . as il- 
lustrated in Figured 

The turbulent nature of an accreti ng protoplanetary disc 
can result in local pr essure maxima (Joh ansen et al.l [2009a; 
iFromang & S tone 200§|). As can be seen from Eq. © these pres- 
sure bumps can locally reduce the headwind the pebbles experi- 
ence. Reductions by 6 A » -0.02 are seen in shearing box sim- 
ulatio ns of the MRI dJohansen et~al . 2009a; Fromang & Stond 
120091) and global simulations dLvra et al.ll2008bl) . We have illus- 
trated the effect of a strong pressure bump, with 6 A « -0.04, in 
Figure [2] 

Since particles face a headwind, they will drift radially and 
azimuthally as 



Tf 



v = - 



T f 2 + 1 



(9) 



(10) 



as s hown by IWeidenschillingl (1 19771) and iNakagawa et al.l 
(1986). The total relative velocity between the particle and the 
core in pure Keplerian rotation is 



AV: 



+ 1 



r? + 1 



-rjVK, 



(11) 
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which is well approximated by Av « r/Vk, since the parti- 
cle sizes we consider, Tf = (0.01,0.1, 1), give us Av/(77Vk) = 
(1.0,1.0,1.1). 

Particles settle in the vertical direction (perpendicular to the 
orbital plane). The particle scale height H p is a balance between 
midplane-directed g ravity and turbulent diffus ion parametrized 
by the coefficient S t dYoudin & Lithwicldl2007h . 



p 
H 



^*0.01, 



(12) 



where we will make the approximation that this holds for the 
pebble size range we consider (Tf = 0.01-1). The turbulence gen- 
erated by the streaming instability (Youdi n & Goodmanll2005l) 
self-regulates the particle midplane density to equal the gas den- 
sity, independent of particle size. Since Z(H p /H)~ l ~ p p /p ~ 1, 
a near-unity midplane density in a protoplanetary disc with 
metallicity Z = 0.01 sets the particle scale height to be H v /H = 
0.01. For turbulence generated through the MRI, a value of 
S t = 0.001 would give a ten times higher particle scale height, 
H v /H « 0.1, for particles of friction time Tf = 0.1. 

The aim of this paper is to investigate accretion onto cores 
of various masses, ranging from the expected initial masses of 
planetesimals to estimated final core masses of gas-giant plan- 
ets. The core mass, or more precisely its gravitational parameter 
GM C , is non-dimensionalized as 



GM C M c (H\ 



MAr) 



-3/4 



(13) 



with G the gravitational constant and M© the stellar mass. Figure 
prelates the dimensionless core mass fi c to the orbital radius in 
the MMSN. Given the core mass, we can assign it an uncom- 
pressed radius of 



R c = 890 



i,2g/cm 



-1/3 



10" 3 M ff 



1/3 



km. 



(14) 



The critical core mass for runa way accretion of a gaseous enve- 
lope is approximately 10 M® (lMizunolll980l) . only weakly de- 
pendent on th e orbital radius outside the terrestrial planet region 
(Rafikov 2006). At 5 AU this m ass corresponds to /i c ^ 1 , as can 
be seen in Figure [3] However, as lHori & Ikomal(l201ll) point out, 
if the envelope can be significantly polluted by heavy elements 
from the accretion of icy bodies, the critical core mass will be 
reduced by up to two orders of magnitude. 

Planetesimals are believed to have i nitial sizes in the 100 - 
1000 km region (I Johansen et al.l 12007k iMorbidelli etai1l2009l) . 
I Johansen et al.l (120121) find the characteristic clump mass by 
streaming instabilities to be yUsi ~ 5 x 10" 6 (see Figure [3]). 

The dynamical equations of the particles, 



dv 

- = -2&k 

at 



? 1 
x v + 3Q^xe x + g c (v - u) , 

tf 



(15) 



are s olved with the Pencil Cod J| in the shearing box approxi ma- 
tion dGoldreich & Tremainelll980t iBrandenburg et al.lll995"l) . A 
Cartesian coordinate system is placed rotating at an arbitrary, 
but fixed orbital distance, with Keplerian frequency Q^. The 
x-axis points radially outwards, the azimuthal direction corre- 
sponds to the v-coordinate and the vertical z-direction is per- 
pendicular to the midplane. The motion of the particles is de- 
scribed by Eq. (ITSb . which includes the acceleration due to the 



1 The Pencil Code can be freely obtained at 
http : //code . google . com/p/pencil-code/ 
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Fig. 3. Contour lines mark the core mass M c (in units of Earth 
masses M©), as function of the dimensionless mass unit yu c = 
(GM C )/(Q^H 3 ) and orbital radius r in AU. The red contour line 
indicates the assumed minimal mass of a gas-giant core. The hor- 
izontal red dashed lines indicate the expected initial seed core 
mass from the streaming instability after planetesimal forma- 
tion by gravitational collapse, //si, and the transition mass fi t (see 
Section 3.2). 



core placed in the centre of the frame and the self-gravity of the 
particles solved for through the Poisson equation, 



V • g c = 4nGp p , 



(16) 



Additionally, it includes the drag force term -jj(y — u), a term 

balancing the linearized gravity and the centrifugal force 3Q\xe x 
and the Coriolis force -2i2 K x v. 

We perform simulations both with and without the backreac- 
tion term of the particles on the gas, with gas backreaction time 
(p v /p)~ l tf. The momentum equation for the fluid elements, 

du 9 9 

— +u-Vu = -2Q K x u + 3Q^xe x - Q^ze z 



VP P P 

+ — —(v-u), 

Pg Pg'f 



(17) 



includes the pressure gradient term, -(l/p p )VP, and vertical 
gravity -Q\z- The continuity equation for the gaseous compo- 
nent of the protoplanetary disc is given by 



!+v.(p«)= 

ot 







(18) 



and we use artificial hyperdiffusivity for the gas to dissipate en- 
ergy on the smallest scales. 

We solve these equations non-dimensionalized by the 
Keplerian frequency £?k, the scale height H of the gas disc and 
po, the gas midplane density. This has the benefit that, when 
interpreting the normalized results, the orbital dependency is 
nearly fully recovered from these parameter's MMSN orbital 
scalings (e.g. Figure [T][2l [3]). The numerical results are however 
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not completely scale-free. When solving for self-gravity with the 
Poisson equation, we must set 



4nGpo 



(19) 



the non-dimensionalised form of the gravity constant G, as an 
initial condition. It shows only a weak dependency on the orbital 
radius, 



r « 0.04 



/ r y/ 4 

\3AU/ ' 



(20) 



in the MMSN. We therefore fix r = 0.04 for the remainder of 
this paper. 

All simulations are performed in a three-dimensional shear- 
ing box, with a fixed particle scale height of H p /H = 0.01. Run 
parameters for all simulations used in this paper can be inspected 
in Table [T] The core is fixed in the origin of the coordinate sys- 
tem, a valid approximation for the range of gas-decoupled seed 
masses we cover. When the escape velocity v esc from the surface 
of the core is small compared the sound speed, variations in the 
gas density can be ignored, as can be seen from the hydrostatic 
equilibrium of an isothermal gas, 



y esc 



d lnp 
dlnr' 



(21) 



We argue this approximation holds up to the largest cores we 
consider in Section 5. All simulations, with the exception of 
le-3_®.l_b (see Table [B, do include the gas drag on the 
particles, but lack the backreaction from the particles of the 
core on the gas. When omitting the backreaction term, the 
gas velocity equals to the sub-Keplerian velocity, u y = -t/vk- 
However, simulation le-3_® . l_b includes the particle's backre- 
ac tion on the gas and follow s the numerical scheme discussed 
in lYoudin & Johans"enl (120071) . When including the backreaction 
term, we also turn on the vertical gravity force for the particles, 
-Q\ze z . 

All runs have sheared periodic boundary conditions in the 
radial direction, but particles crossing azimuthal boundaries 
get removed from the simulation domain (with exception of 
le-3_®. l_b), in order to avoid accretion of particles already 
focused from their first passage past the core. We have also 
run simulations including collisions, with the scheme discussed 
in iJohansen et al.l (120121) . and found no measurable difference 
in the accretion rates on the seed core masses. We therefore 
omit collisions from the simulations in this paper. Implications 
and limitations of the simulation set up are further discussed in 
Section 5. 



3. Results 

Inspection of the particle's momentum equation, Eq. ([T5l) , re- 
veals an important length scale. The Hill radius, 



rH 



GM C 



\3Q 2 



= ( 



,1/3 



(22) 



is set by the gravitational competition between the acceleration 
towards the core and the stellar tidal field in the radial direc- 
tion. At a separation rn from the core, the orbital time around 
the core approximately equals the orbital time around the star, 
2nQ^}. The Hill sphere's radius grows linearly with the orbital 
radius rn oc r, placing more material in the gravitational region 
of influence of the core. 



3.1. Drift accretion 

When ignoring the stellar tidal field and the Coriolis force, the 
Bondi radiu£| 



GM C 
Av 2 



(23) 



marks the outer point at which particles approaching the core 
with relative velocity Av get significantly gra vitationally de- 
flected (> 1 rad, e.g. iBinney & Tremainel 119871) . Ignoring the 
stellar tidal field is a valid approximation before the core mass 
grows to the point where the Bondi radius becomes comparable 
to the Hill radius (see Section 3.2), and we can associate a core 
mass, 



1 Av 3 



(24) 



with this transition. Masses with M c < M t are in the drift regime 
and pebbles embedded in the gaseous disc approach the core 
with a mean velocity comparable to the gaseous headwind the 
core experiences, Av « tjvk- 

Depending on the balance between the gravitational attrac- 
tion of the core and the drag force the pebbles experience, a par- 
ticle can be pulled from the mean gas flow and accreted if it dis- 
sipates enough energy while being deflected. When the gas-free 
core-crossing time associated with the Bondi radius, 



t - rB 
Av 



(25) 



is similar to the friction time tf, the drag force will cause all 
pebbles within the Bondi radius to spiral inwards. However the 
effective accretion radius, the drift radius r<j, shrinks with respect 
to the Bondi radius, when tB ~ tf is not satisfied. 

When tB > tf, the particle under consideration is strongly 
coupled to the gas. In this limit, only grazing particles deflected 
on time- scales shorter than the friction time get pulled out of the 
flow. If we let g denote the gravitational attraction due to the 
core's mass, the condition 



Av 

U = — <tf 



(26) 



needs to be satisfied for accretion to occur. Since the deflection 
time t g is given by Avr 2 /(GM C ) = (r/re) 2 tB, the effective drift 
accretion radius is given by 



t B 

r d = - 



-1/2 



(27) 



in the strong c oupling limit. Th i s radiu s corresponds to the set- 
tling radius in Ormel & Klahr (2010) and is also equal to the 
radius found by Perets & Murray-Clayl(l201ll) where drag forces 
shear apart bound binaries in the Epstein regime. We verified this 
power law by numerically integrating orbits of test particles in 
the 2-body problem including drag, 



dv x /Av 

dt/t B 
dvy/Av 

dt/t] 



B 



_/^B\ 3 X_ _ tB Vx 

\ r / r B tf Av 9 

I r / r R tf \Av / 



(28) 
(29) 



2 Note that we define the Bondi radius with the square of the relative 
velocity between core and particle in the denominator, and not the sound 
speed squared, which is also found in the literature. 
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Table 1. Characterizing parameters of all simulations used in this paper. All simulations are performed with 128 3 grid cells resolution 
in stratified shearing boxes, with particle scale height H p /H - 0.01 and metallicity Z = 0.01. The first column gives the name of the 
simulation, followed by the characterizing parameters: the core mass yu c , particle size Tf , headwind parameter A and side length L of 
the cubic simulation domain. The last column indicates whether the simulation includes the backreaction (BR) term of the particles 
on the gas, or not. 



name 




Tf 


A 


L/H 


BR 


le-6_0.01 / le-6_0.1 


lO" 6 


0.01/0.1 


0.05 


0.01 


No 


le-6_0.1_0.03 


10" 6 


0.1 


0.03 


0.01 


No 


2.5e-6_0.1_0.03 / 2.5e-6_0.1 / 2 . 5e-6_0. 1_0.07 


2.5 x 10" 6 


0.1 


0.03 / 0.05 / 0.07 


0.02 


No 


le-5_0.01 / le-5_0.1 / le-5_1.0 


10" 5 


0.01/0.1/1 


0.05 


0.04 


No 


le-5_0.1_0.03 / le-5_0.1_0.07 


10" 5 


0.01/0.1/1 


0.03/0.07 


0.04/0.01 


No 


le-4_0.01 / le-4_0.1 / le-4_1.0 


10" 4 


0.01/0.1/1 


0.05 


0.128 


No 


le-3_0.01 / le-3_0.1 / le-3_1.0 


lO" 3 


0.01/0.1/1 


0.05 


0.32 


No 


le-3_0.1_b 


10" 3 


0.1 


0.05 


0.2 


Yes 


le-2_0.01 / le-2_0.1 / le-2_1.0 


io- 2 


0.01/0.1/1 


0.05 


0.64 


No 


le-l_0.01 / le-l_0.1 / le-l_1.0 


10" 1 


0.01/0.1/1 


0.05 


1.28 


No 



where we non-dimensionalized the particle equation of motion, 
ignoring disc dynamics. This is a valid approximation in the drift 
regime, where <?c Q^} . Sample orbits can be investigated in 
the inset of Figure ffl which shows the maximal particle-core 
separation leading to capture. The drift radius for strongly cou- 
pled particles falls of as oc (7b Af)~ 1/2 as predicted. Particles with 
tf ~ tB, get efficiently accreted within a Bondi radius from the 
core. 

Particles weakly coupled to the gas with respect to low-mass 
cores (7b < tf) are less aided by drag as they get deflected by 
the core. As seen in Figure [4] a rapid fall-off occurs for particles 
with tf « 10 2 ^b- The orbits in the inset show these particles to 
be gravitationally scattered, similar to the case were no gas drag 
is present. Here, the physical radius of the core becomes rele- 
vant, since accretion now occurs through gravitational focusing 
of particles on the core's surface, which we have not taken into 
account in Figure @1 

The accretion rate in the drift accretion regime is given by 



M d = 7rp p r d Av, 



(30) 



when rd is smaller than the particle scale height H p . A represen- 
tative simulation in this regime, performed with \i c - 10" 5 , is 
illustrated in Figure [5] Pebbles drift with a sub-Keplerian veloc- 
ity past the core and those entering the Bondi radius, here well 
inside the Hill radius, feed the growth of the embryo. Note that 
when rd ~ rb, the core growth scales faster than exponential with 
mass, as Md oc M 2 . Figure [6] shows the accretion rates calcu- 
lated from simulations le-6_® . 1 , 2 . 5e-6_® . 1 and le-5_0 . 1. 
Particles with friction time Tf = 0.1 closely follow the maximal 
drift accretion efficiency, 



/i d 1 P P / //c 



(31) 



with rd ~ tb. However the low-mass core in run le- 6-0.1 
comes close to the weak coupling limit and sees its accretion 
rate reduced. 

We can envisage two effects reducing the accretion rate, if we 
were to continue to ignore the stellar tidal field even for higher 
mass cores. Firstly, when the core enters the strong coupling 
limit, growth slows down to exponential, Md oc t~ l r 2 oc M c . 
Secondly, when the accretion radius becomes comparable to the 
particle scale height, the appropriate expression for the accretion 
rate is given by 



M d = 2r d 2 p Av, 



(32) 



10° 
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Fig. 4. Accretion efficiency in the weak and strong coupling 
regime. When the Bondi time t^ = GM/Av 3 is equal to the 
particle's friction time tf, the drift accretion radius r^ peaks and 
equals the Bondi radius r%. For a particle of fixed size, the ratio 
tB/tf on the horizontal axis increases as the core mass grows in 
time. When particles are strongly coupled to the gas (fe > tf), 
with respect to the gravitational attraction of the core, the drift 
radius decreases as r$ ~ (W4)~ 1/2 (the analytical scaling of 
Eq. (l27l) is indicated with a full grey line). Near t^/tf « 10" 2 the 
drift radius rapidly decreases. The inset shows particle trajecto- 
ries (grey curves) in this regime, which can be compared with 
those at fe/tf - 1 (black curves). Where the former are simply 
gravitationally deflected, in the latter case we see that particles 
inside the Bondi radius (marked by a red circle) are accreted by 
the central point source. 



where X p is the particle column density. When r& 
exponential growth Md/M c = 2E p /Av. 



tb, we get 
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0.02 




-0.01 



0.00 

x/H 



0.02 



Fig. 5. Accretion of pebbles with Tf = 0.1 by the central core (yu c = 10" 5 ) in the drift regime. The color coding shows the local 
particle surface density £ p , normalized by the average particle density (^ P ), in the simulated shearing box with Z = H p /H = 0.01. 
Marked as a white dot, the central seed core can be seen in the first panel. Both the drift and Hill radii are plotted as white circles. 
The second panel illustrates the drift of the particles (A = 0.05) and creation of an accreting particle wake. In the third panel, 
particles within the drift radius r& « r% are accreted. Particles further out may be carried out of the box by the sub-Keplerian gas or 
by the Keplerian shear. In the final panel the headwind has blown most pebbles past the core, with only a minority accreted. 



3.2. Hill accretion 

When the core mass grows to the point where the Bondi radius 
r B oc Ml is comparable to its Hill radius rn °c M* /3 (or identi- 
cally vh = Av or tn/tf = r^ 1 ), it will cross the transition mass, 



1 Av 3 



3GQ K 



3x 10" 



A \ / r 
O05 IsAij) 



(33) 



defined earlier and see a change in pebble accretion mechanism. 
The dimensionless form of the transition mass, 



// t = 7x 10" 



0.05 



(34) 



scales as the cube of the headwind parameter A (see Figure©. 
The Hill radius now sets the maximal impact parameter from 
which particles can be accreted. When M c > M t , pebbles at the 
edge of the Hill sphere approach the core with relative velocity 



(35) 



The Keplerian shear v = -(3/2)Q^x domi nates o ver the head- 
wind in the Hill branch, since vh/Av = V^b/^h- The inverse 
Keplerian frequency, -Q" 1 , is the gravitational crossing time- 
scale at the Hill radius, independent of core mass. For particles 
with friction times close to the orbital time-scale (Tf = 0.1-1), all 
particles entering the Hill sphere will be accrete, as illustrated in 
Figure H Here we present particle orbits obtained from the Hill 
equations including drag 



dv x /vn 
dt/Q^ 1 

dvy/vu 



+2— + 3 31 — ) 

vh rn \r H / r H r f v H 



v H 



= -2— -3 — —--{— + - — 



y 



l 



3 x 



r H T f \v H 2r H 



(36) 
(37) 



where r - ^x 2 + y 2 is the particle-core distance. Similar to the 
drift case (Section 3.1), when the gravitational deflection time 
(here independent of the core mass and ~ Q~ l ) is similar to tf, 
enough energy will be dissipated during the approach to mediate 
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Fig. 6. Accretion rate fi/po as function of the initial core mass //o. 
Theoretical curves for the branches corresponding to drift and 
Hill accretion are plotted as respectively dashed and full lines 
in grey. The grey full circle marks the transition mass. Black 
crosses represent the simulated results in a stratified shearing 
sheet, with A = 0.05. Triangles correspond to simulations with 
modified A, a grey upwards pointing triangle corresponds to A = 
0.07 and a downwards pointing triangle corresponds to A = 0.03. 
In the Hill branch, the position of both triangles lie on top of the 
black crosses and are omitted for clarity. The grey square shows 
the result of simulation le-3_® . l_b, which includes the particle 
backreaction. 



the accretion of the pebbles within the Hill sphere. The accretion 
rate is then given by 

M H = 2r H 2pVH K M c /3 ' ( 38 ) 

since for the core masses under consideration, M c > M t , 
we accrete the total particle surface density, S p (ru > 
H v ). This growth mode is confirmed by our numerical sim- 
ulations le-4_0.1, le-3_®.l, le-2_®.l, le-l_®.l and 
le-3_1.0, le-2_1.0, le-l_1.0. They can be inspected in 
Figure [6] and [8] An example of a simulation with a seed mass 
accreting at the Hill rate is illustrated in Figure Here, particles 
of Tf = 1 entering the Hill sphere drive the growth of the core 
of mass p c = 10" 2 . Accretion occurs through a particle disc, 
as was previously resolved in high-resolution 2D simulations by 
I Johansen & Lacerdal (l2010h . 

Note that in the classical scenario of planetesimal accretion, 
one never captures objects from the full Hill sphere, but only 
by a fraction a l/2 rn, with a « r c /rn s et by the physi cal ra- 
dius of the core r c . In the terminology of Rafikovl (1201 ll) . slow 
accretion of planetesimals between the shear- and dispersion- 
dominated dynamical re gime, from a part of the parti cle scale 
height H v = v/Q^ (e.g. Do dson-Robinson et al.ll2009b . goes as 
approximately M « nar^p^v « arn^pVu ~ aM\\. At 5 AU, this 
gives a reduction in the accr etion rate of a g 10~ 3 (r/5 AU)" 1 , 
for a standard solid density (iGoldreich et al.l 12004). Accretion 
of planetesimal fragments from a thin midplane, as discussed 




-2 



-1 







Fig. 7. Trajectories for particles with dimensionless friction time 
Tf = 0.01, 0.1, 1 obtained from the 2D Hill equations including 
gas drag. Pebbles with Tf = 0.1 and impact parameters below 
a Hill radius efficiently get accreted. Larger particles of Tf = 1 
are pulled in from wider separations, but cores lose particles on 
horseshoe orbits. Particles strongly coupled to the gas, with Tf = 
0.01, need close encounters well within the Hill sphere in order 
to fall onto the core. 



m Rafikov (2004) is more efficient, with the accretion rate be- 
ing proportional to V^Mh. At 5 AU, this limits the growth by 
^[a « 3 x 10" 2 (r/5 AU)" 1/2 . Thus accretion of pebbles at a rate 
Mu from the entire Hill sphere is extremely efficient, compared 
to the classical gas-free case. 

Only the smallest particles we consider, with Tf = 0.01, have 
an accretion efficiency that is less than optimal in the Hill branch 
(Figure [8]). Similar to the case of the strongly coupled particles 
in the drift regime, accretion requires the gravitational deflection 
time to be shorter than the friction time, as previously expressed 
in Eq. (l26l) . The relative velocity for particles approaching the 
Hill sphere is set by the Keplerian speed Av « Qr. This allows 
us to rewrite the accretion criterion as 



Qr < tf, 

GM 

which gives us an effective accretion radius 

. 1/3 



(39) 



(40) 



In this regime, Mn, e ff 



OC T 



2/3 



which compared to particles of 



Tf = 0.1 would give a reduction of the accretion rate by 
as can be seen from comparing Figure and [8] 



0.2 



3.3. Influence of headwind reduction and particle feedback 

In the above discussion, we have kept the relative velocity be- 
tween core and the gas disc constant at A = 0.05. As pre- 
viously mentioned, we have ignored the presence of pressure 
bumps, local extrema in the radial pressure force resulting in re- 
gions of a reduced headwind, as well as extreme orbital distances 
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Fig. 8. Normalized accretion rates, /* C M), for different particle 
sizes with friction time Tf = 1 and Tf = 0.01, as indicated by 
respectively asterisks and diamonds. The grey curves correspond 
to those shown in Figure [6] and similarly the transition mass is 
indicated by a full grey circle. In the Hill branch, larger particles, 
Tf = 1 , get accreted as efficiently as particles with friction time 
Tf = 0.1 (see Figure©, but in the drift branch they never get 
accreted at the full drift rate. On the other hand, small particles, 
Tf = 0.01, get efficiently accreted in the drift branch, but less so 
in the Hill regime. 



where A can change significantly. While the accretion rate Mr 
in the Hill regime is insensitive to A, the Bondi branch up to 
the transition mass M t is not. Figure [6] illustrates the effect on 
changes in A for various core masses (A = 0.05 ± 0.02). For a 
core mass accreting approximately from the full Bondi branch, 
Eq. (OTI) indicates that the accretion rate will be modified by a 
factor (0.05/0.03) 3 « 5 for A = 0.03 compared to accretion with 
A = 0.05. This is in agreement with the measured accretion rates 
for ji = 2.5 x 10" 6 (2 . 5e-6_® . 1_® . 03). The increase is reduced 
for le-5_® . 1_® .03 where r B grows to Ru, and increased for 
le-6_®.l_®.®3 where t^/tf grows sufficiently out of the weak 
coupling limit. Overall we see that even weak pressure bumps 
decreasing A by 0.02 lead to more rapid accretion. Vice versa, 
for A = 0.07 we expect a reduction of the accretion rate by a 
factor (0.05/0.07) 3 « 0.4 and we measure similar, but lower ac- 
cretionrates (le-5_0. 1_0.03 , 2 . 5e-6_® . 1_0 . 03). 

In simulation le-3_® . l_b, we have departed from a smooth 
gas velocity profile, by including friction on the gas and letting 
turbulence develop by the streaming instability. After approxi- 
mately 20 orbits, we place the seed core mass in the center of 
the simulated domain (Figure [TO]). The measured accretion rate 
does not deviate measurably from the case not including particle 
backreaction, as can be seen in Figure This indicates that our 
results are robust for the Hill branch, even in a turbulent envi- 
ronment. More simulations should be carried out in the future 
to verify the validity in the drift regime. However, this requires 
very high resolution in order to resolve both the Bondi radius 
and the streaming instability wavelength simultaneously. 



4. Implications for gas giant growth 

Having numerically confirmed the pebble accretion rates for low 
to high core masses, we can extrapolate our results and find the 
time necessary to grow a core to the critical mass needed to at- 
tract its gaseous envelope. 

In the drift regime, the accretion rate implies a growth time- 
scale of 



rM c 
Jmq 



M d l dm 



Av 3 
np p G 2 



M, 



-l 



(41) 



to reach a core mass M c from an initial seed mass Mo <?c M t . 
This lower limit, since we assume optimal accretion from the 
full Bondi radius, sets the time until the accretion rate blows up 
hyperbolically as M c oc (Af<j - t)~ l and thus does not depend on 
the final mass we wish to reach. However, from our numerical 
results we know that this growth is not sustained and turns off 
to the Hill branch. The characteristic time to reach the end point 
for drift accretion, the transition mass M t , is given by 



At d * 8 x 10 e 



A 
O05 



Pp/PO 

0.01 



M 



10" 5 ]VL 



f r y 



yr, (42) 



which at 5 AU is comparable to the gas disc lifetime. If parti- 
cles sediment to the midplane, the ratio of the particle to the gas 
density, p p /po, would be of order unity. However, past 0.5 AU, a 
small seed core mass, Mo = 10" 5 M®, accretes non-sedimented 
particles of size Tf < 0.01 most efficiently. For reasonable values 
of the local headwind, A = 0.03 - 0.07, growth is too slow to 
form cores large enough to enter the Hill accretion regime, and 
without pressure bumps A is even larger at wide stellar separa- 
tions. 

Hill accretion on the other hand, has a growth time- scale of 



rM cr 

At n = 



M H dm 



3 5 /W /3 
- M 1/3 

crit ' 



2G 2 ' 3 l, n 



(43) 



which is only weakly dependent on the critical mass for gas en- 
velope attraction M cr i t and independent of the transition mass, 
when M t <^c M C n t . The core growth when accreting pebbles is 
fast in this regime, at 5 AU the critical mass is reached after 



A^h ~ 4 x 10 4 



M C1 



10 M e 



1/3 



Gau) 



yr. 



(44) 



Furthermore, the growth time-scale A^h scales linearly with or- 
bital distance r, as opposed to quadratic in the drift regime. This 
makes core formation possible in distant regions of the proto- 
planetary disc. Also, note that Hill accretion rate is maintained 
for a single particle size with friction times Tf ~ 0.1-1, indepen- 
dent of the core mass, as opposed to the Bondi regime where one 
unrealistically needs to maintain « tf to maintain the maximal 
accretion rate. 

Figure [TT] shows the core growth in both regimes and the 
dependency on the orbital distance. We conclude that fast core 
growth is possible through pebble accretion, provided that the 
initial seed mass for the core is above the local transition mass. A 
sufficiently large embryo can only be grown by drift accretion in 
pressure bumps with low headwind, A < 0.05, or be the result of 
planetesimal formation by gravitational collapse after concentra- 
tion by e.g. streaming instabilities (see also Figure[3]). It is inter- 
esting to note that both Ceres and Pluto have less than the critical 
mass needed for fast Hill accretion, which might explain why 
they failed to grow to gas or ice giants. Indeed, one can make 
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Fig. 9. When the core is massive enough, it can efficiently accrete particles with Tf = 1 entering its Hill sphere (indicated by the 
white circle in the first panel). On this scale, the Keplerian shear dominates over the relative velocity difference between the gas and 
Keplerian velocity (A = 0.05). Accretion seems to occur through a particle disc, visible after a steady state has set in (panel 2 to 3). 
In the last panel, accretion and Keplerian shear have removed most particles in the box. The color coding is similar to Figure[5] The 
simulation was performed in a stratified shearing sheet box with Z = H p /H = 0.01. 



the assumption that only those planetesimals that formed early 
enough in the high-mass tail of the initial planetesimal mass dis- 
tribution could serve as the seed for gas-giant cores. 

As an illustration of the rapid core growth by pebble accre- 
tion in the Hill regime, we compare it to the core growth time 
for planetesimal accretion in Figure [12] As discussed in Section 
3.2, the inability to accrete solids from the entire Hill sphere, as 
opposed to pebble accretion, leads to significant longer core for- 
mation times, in conflict with the observed dissipation time of 
protoplanetary discs. 

5. Discussion 

We discuss here the assumptions and limitations of our results. 

Midplane layer thickness. One component of the pebble ac- 
cretion scenario is the presence of a thin particle disc (H v = 
0.01//). This low particle scale height is expected from tur- 
bulence driven by streaming instabilities, independent of parti- 
cle size, as discussed in Section 2. A moderately higher par- 
ticle scale height, as may be the case for turbulence caused 



by the magnetorotational instability, can result in a situation 
where th < H p past the transition core mass. This would re- 
sult in a temporarily reduced accretion rate, by a factor rn/H p = 
(H p /H)- l (r n /H) = (l/3) 1/3 (# p /#)-y /3 , until the Hill radius 
grows beyond the particle scale height. 

Particle size. The assumption of a single particle size in our 
simulations can be criticised, but as discussed in the introduc- 
tion, observations of protoplanetary discs allow a large fraction 
of the solid mass t o reside in the pa rticle size range that we con- 
sider, Tf = 0.01-1 dWilner et al .120051) . A large abundance of par- 
ticle s larger than pebbles i s not expected from c oagulation mod- 
els (iBlum&Wurml 120081: | Braueretal . 2008; Windmark et al. 
12012^ However, as particles approach the core their icy com- 
ponent might sublimate; as friction would heat the particles, es- 
pecially when a denser envelope starts forming around the core. 
It would be interesting to take this size-diminishing effect into 
account in a further investigation. On the other hand, particles 
might grow larger. In higher metallicity environments streaming 
instabilities become so effective in clumping solid material that 
one can fear particles to grow past the pebble size. However, 
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Fig. 10. Including the backreaction friction force on the gas, a core of [i = 10" 3 accretes particles of Tf = 0.1 (third and last 
panel), even in presence of turbulence caused by the streaming instability. The turbulence has first been given 20 orbits (-126 Q~ x ) 
to saturate (first and second panel), before the core is inserted. From the full Hill sphere, r^/H « 0.07, pebbles are attracted in a 
prograde motion to the core (third panel). As before, the color bar in the first panel gives the particle surface density. The insets 
show the azimuthally averaged particle density and vertical extent of the particle layer. For clarity, in the insets the color coding 
covers a twice as wide range in particle overdensity, compared to the surface density plots. 



we do not see this particle clumping in our simulations includ- 
ing the gas drag backreaction at the metallicity we consider 
(Z = 01). Strong c lumping requires Z > 0.02 (iJohansen et al.l 

|200^Ib^^ 

Gas structure. For the lower seed masses discussed in the 
paper, we previously argued (Section 2) that the gas density 
changes around the core are small. In the Hill regime, the ra- 
tio vljcl « 2.3 x 10y /3 (r/AU) (in the MMSN for standard 
solid density) can exceed unity for the highest core masses and 
the effects of an envelope should be ta ken into account. But, as 
also argued by lOrmel & Klahr I (l2010h . even if the direction of 
the flow moderately changes on scales within the Bondi radius 
due to stratification near the core, only particles with tf <*c tj$ 
could be affected by it. Since these particles are too strongly 
coupled to the gas for accretion to take place in the first place 
(strong coupling limit), ignoring the core's feedback on the gas 
is justified. 

Keplerian orbits. In our analysis we assumed the core to 
be on a circular Keplerian orbit. The relative velocity between 



the core and the gas in Keplerian rotation could be significantly 
modified if competing cores would get excited by repeated close 
passages. However, as opposed to classical planetesimal growth, 
in our scenario gas damps the small particles and dynamical 
friction prevents the excitation of larger bodies, similar to the 
oligarchic growth regime. We do ignore gas-driven type-I mi- 
gration of the core , important for core masses over 0.1 M© 
(iTanaka et alJl2002h . 

Random particle speed. In our simulations, particles ap- 
proach the core in equilibrium with the gas flow. Particle in- 
teractions with the core last at most of the order I?" 1 , as in 
the Hill regime. Small particles (Tf < 1) are coupled to the 
gas on similar time-scales. The passage of the core is quickly 
erased for the non-accreted particles, even when ignoring ra- 
dial drift and turbulent diffusion. The core only catches up with 
the deflected particles after approximately £ pa ss = 27rrAv _1 = 
(27r/Q K )(r/H)A- 1 « 10 4 (r/AU) _1/2 Q^. However, in the Hill 
accretion regime all particles that can be deflected are accreted, 
and it is the radial drift and diffusion of particles that fill up the 
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Fig. 11. Core growth as function of time, plotted for various or- 
bital distances (0.5, 5 and 50 AU). The drift branch, marked by 
grey solid lines, assumes an initial core mass of Mo = 10" 5 M® 
and A = 0.05. The drift growth continues until the transition 
mass Af t is reached (marked by a full grey dot). Accretion con- 
tinues through the more efficient Hill branch, drawn in black. 
For clarity, we start the Hill growth from the transition mass at 
time t = Oyr, instead of continuing from the time where drift 
accretion comes to a halt. The masses of Ceres and Pluto (lo- 
cated at respectively 2.7 and 39 AU) are marked on the vertical 
axis for reference. The grey dotted curves correspond to clas- 
sical planetesimal accretion (PA), where the faster gro wth cor- 
respo nds to 2D accretion of planetesimal fragments (Rafikov 
200 4() and the slower to 3D accretion of planetesimals (e.g. 
Dodson-Robinson et al.l l2Q09h . Note that drift accretion time- 
scale at 50 AU takes more than 10 8 yr and its transition mass 
point is not plotted. 



feeding zone. Diffusion can be rapid, since the diffusion time 
associated with closing the Hill sphere t\\Q ~ R^/(6 t H 2 ) ~ 
(l/3) 2/3 £ t -1 yU 2/3 , is of order unity for a protoplanetary disc with 
£ t = 0.01. 

Particle drift. When the drag force responsible for radial 
drift is too small, particles could get trapped in mean mo- 
tion resonances with the core. Weidenschilling & D avisl (Il985l) 
studied large, Tf > 1, particles in the Stokes drag regime, 
and argued that particles smaller than these sizes feel large 
enough drag forces to escape resona nt trapping around a J upiter- 
mass planet at 5 AU. As shown by iTanaka & Idal (1 19971) inclu- 
sion of mutual planetesimal interactions breaks down the reso- 
nances, but dust gap formation still occurs for large planetes- 
imals, where gas drag changes the semi-major axis of the the 
planetesimals after scattering with the protoplanet. The maxi- 
mal particle size unaffected by particle trapping seems approxi- 
mately inversely proportional to th e planet's mass, whic h is also 
seen in simulations performed b y Haardekooper (2007). In fact, 
Weidenschilling & Davis (1985)' argue that small pebbles are the 
only size that can be accreted by the core, since trapped larger 
planetesimals get dynamically excited and will be ground to 
fragments, which in their turn are capable of escaping the res- 
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Fig. 12. Time needed for core growth up to 10 M© at various lo- 
cations in the disc. The solid black line gives the formation time 
of the core for pebble accretion in the Hill regime, while grey 
lines give the time needed to form the critical 10-Earth-mass 
core by planetesimal accretion. The dashed grey line represents 
planetesim al fragm e nt acc retion from a thin midplane layer, as 
studied by Rafikovl (120041) . The red shaded area shows the ap- 
proximate time interval in which the protoplanetary disc loses 
its gaseous component and enco mpasses for example th e esti- 
mated age of gas giant LkCa 15b ( Kra us & Irela nd 2012). Core 
formation needs to occur before this time. 



onance. This picture is confirmed in simulations performed by 
iLevison et al1(l201Ql) . 

Dust gaps can open up before th e core is massive enough 
to create a gap in the ga s disc itself (iPaardekooper & Mellema 
120061) . iMuto & Inutsukal d2009h analytically show that the core 
has to be over a critical mass, 



(H \ 
— J ~ 1, 



(45) 



for particles of Tf < 1 in order for a dust gap to emerge. Past r sb 
1 AU, ii sb 1 is consistently above 10 M®, the critical core mass 
for gas and ice giants (see Figure [3}. Particles thus always drift 
radially fast enough to replenish the feeding zone of the core. 
Indeed, if the drift rate is set by R& = 2^rAvE p , the requirement 
Rd > Mu recovers th e above criterion, Eq . (H5V At the same 
time, as pointed out by lOrmel & Kobayashil (120121) . the particle 
drift can also be responsible of clearing up the entire reservoir of 
available pebbles in the disc. 

Terrestrial planet formation. Growth at small orbital dis- 
tances, r < 5 AU, is remarkably rapid in the pebble accretion 
model. Formation of rocky planets and possibly in situ forma- 
tion of gas-giant planets in the terrestrial planet region seems 
problem-free from the perspective of the accretion rate. The 
growth time- scales for both the drift and Hill accretion branch 
shrink to approximately 10 5 yr at Earth-like separations from 
the host star. This could indicate that even terrestrial planet for- 
mation occurs rapidly during the gaseous disc phase. However, 
closer to the star the amount of material in an annulus of Hill- 
radius-width is small and the isolation mass by gap formation is 
lower. Also the optimally accreted particle size is large, around 
10 cm, and ices are not available. 



6. Conclusions 

In this paper we have demonstrated that accretion of pebbles 
makes rapid formation of gas-giant cores possible. The growth 
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time-scale to reach the critical core mass for gas accretion is 
reduced by three orders of magnitude at 5 AU and four orders 
of magnitude at 50 AU, compared to the planetesimal accretion 
rate in between the shear- and dispersion-dominated dynamical 
regime. Compared to accretion of planetesimal fragments from 
a thin layer the formation time is shortened by approximately a 
factor 30 at 5 AU and a factor 100 at 50 AU. This is further sup- 
port for the core accretion scenario, because cores can form by 
pebble accretion before gas dissipation after 1-10 million years, 
even at large orbital radii. 

We can summarise the main numerical results as follows. 
Our simulations show gas drag to be a necessary ingredient for 
fast pebble accretion by the growing core. Omnipresent pebbles, 
particles with friction time around Tf « 0.1, are ideally suited 
for core growth. They are weakly enough bound to the gas to 
feel the gravitational pull from the core, but strongly enough to 
deposit their kinetic energy through drag forces, when passing 
the core. Low-mass cores, cores below the transition mass cor- 
responding to a body of radius larger than approximately 1 , 000 
km, can accrete small particles drifting with the sub-Keplerian 
gas velocity past the core, but this process is slow, even in pres- 
sure bumps with reduced headwind and particle settling in a thin 
mid-plane layer. However growth in this regime could be im- 
portant for a seed planetesimal formed just below the transition 
mass, where the accretion rate is high. Higher-mass cores can 
e fficiently attract pebble s from the full Hill sphere, as was fo und 
bv I Johansen & Lacerdal (l2010h and lOrmel & Klahr I (l2010h . In 
this regime, the optimally accreted particle size is independent 
of the core mass, in contrast to the drift regime, where the par- 
ticle size with the highest accretion rate increases linearly with 
radius of the growing core. 

For the pebble accretion mechanism to be rapid, a significant 
fraction of the solid density needs to be in the form of pebbles 
close to the midplane, and some planetesimals need to form with 
sizes of 1000 km or larger. Theoretical models of planet forma- 
tion show that these large seeds of approximately Ceres-size can 
form by self-gravity after c lumping by the streaming instability 
(I Johansen et al.l2007[|2012l) . so the pebble accretion scenario fits 
well with the formation of planetesimals by self-gravity. 

The conditions for fast core growth are supported by ob- 
servations. A large reservoir of pebbles is i nferred in observa- 
tions of many young protoplanetary discs (IWilner et al.l 120051: 
iRodmann et al] [2006). Additionally, studies of the collisional 
evolution of the asteroid belt show that la rge asteroids must hav e 
formed early when gas was still present (Morbi ddli et al.ll2009l) . 
The early disappearance of mm-d ust in protoplan etary discs on 
time-scales shorter than 1 Myr (iLee et al.l 12011 1) can be con- 
tributed to fast particle growth and rapid core formation. In fact , 
we see that gas giants for m both rapidly (Kraus & Irelandll2(5l2l) 
and at large orbital radii dMarois et al.ll2010l) . Our results predict 
that ice and gas giant planets, detectable with direct imaging sur- 
veys, will be abundant around young (~1 Myr) stars. 

Further exploration of the pebble accretion mechanism by 
simulations with particle backreaction, around higher or lower 
metallicity discs in larger simulation domains, are needed to 
show the robustness of rapid core growth. Preferably, global 
simulations should be developed, including gap formation, ra- 
dial drift, mean motion resonances and multiple cores, in order 
to get a full overview on the implications on fast core growth. 
When the core mass starts to approach the critical 10 Earth- 
masses for envelope attraction, we see a prograde particle disc 
emerge. Studying these discs, possibly the birthplaces of the reg- 
ular satellites of gas giants, in conjunction with the emergence of 



a dense envelope around the core, will teach us more about the 
early growth and gas accretion thermodynamics of gas giants. 

Pebble accretion provides us with a viable pathway to rapid 
formation of gas-giant and ice-giant cores. The necessary first 
step of the core accretion scenario can occur even at wide stellar 
separations, well within the lifetime of gaseous protoplanetary 
discs. 
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